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Section 1 
ALGEBRA 


1.1. POWERS OR INDICES 
If n is a positive integer, 

x” =X XX XX... xX xX acontinued product of n factors 
Otherwise, provided x + 0, 

(yn =0,x° = 


(b) n = —p, where p is real and positive, x-? = = 


(c) n = p/q, where p, q are both integers, x?!* = %/x? 
If x = 0, 0" = 0 for all values of n. 


Laws of Indices 


x™ XX xX x mtn x™ > x? = m0 (x*)* = xm 
1.2, LOGARITHMS 


Derinition: If x = a°, then log, x = c. 


log (xy) = log x + log y 

log (x/y) = log x — log y 

log (x") = n log x 

log (*\/x) = (1/n) log x 

logiox = ae X loge x = log, x/log. 10 = 0:4343 log. x 
(e = 2-718. . .) 


1.3. IMPORTANT FACTORS 


x? — a? = (x + a) (x — a) 
x? + a® = (x + a) (x? F ax + a?) 
x? + y? + 2° — 3xyz 
= (x + yt z)(X? + y? + 2? — xy — yz — zx) 


1.4. QUADRATIC EQUATIONS 

If the quadratic equation ax? + bx + c = 0 has roots «, 8, then 
x? + (b/a) x + cla = (x — a)(x — B) = x? — (a + B)x + of 
Sum of roots = « + 8B = —b/a 
Product of roots = «8 = c/a 
a, B = [—b + (6? — 4ac)]/2a 

Discriminant D = b? — 4ac 


D > 0, roots real and unequal 
D = 0, roots real and equal 
D < 0, roots complex and unequal (complex quantities, p. 23) 


1.5. REMAINDER THEOREM 

If n is a positive integer and ap, Gn-1,. . ., 41, Ao are Constants, 
SF) = ayx" + Qn-1 xn-i +... t+ax+a 

(x) is a polynomial of degree n. 
If f(x) is divided by (x — c), the remainder is f(c). 
Thus (x — c) is a factor of f(x) if f(c) = 0. 

1.6. RATIO AND PROPORTION 


The ratio of a to b (written a:b) is denoted by the fraction a/b. 
a, b, c and d are in proportion if a/b = c/d 
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If a/b = c/d then 


a+b cod 
ad = be co 
a+b ctd la+mb_Ic+md 
a—b c—d pa+qb  pa+qd 


(where /, m, p, q are real numbers). 


1.7. VARIATION 


If k is a constant, 

Direct variation: if y o x, then y = kx 

Inverse variation: if y o 1/x, then y = k/x 

If z oc x when y is constant, and if z oc y when x is constant 
then z « xy, ie. z = kxy when both x and y vary. 


1.8. SIMPLE SERIES 


An Arithmetic Progression is a series of numbers such that any 
term is derived from the previous term by adding to it a constant 
quantity, called the common difference (d). 
Let a be the first term of the series, then the nth term is 
l=a+(n—1)d 
and the sum of 1 terms is given by 


S, = a+ (a+d)+(a+2d)+?..+(@+n-— 1d) 
= $nQQa +n — 1d) = n(a +I) 
A Geometric Progression is a series of numbers such that any 
term is derived from the previous term by multiplying the latter 


by a constant quantity called the common ratio (r). The nth 
term is ar"~* where a is the first term. 


Sx, =at+ar+ar’+...t¢a° =al—r/i—-n 
The sum to infinity = So =a/(l—r) (|r| < 1) 


Harmonic Progression. A series of numbers a, b, c,. . . will be 
in Harmonic Progression if their inverses 1/a, 1/b, 1/c, . . . are 
in Arithmetic Progression so that 


l/c — 1/6 = 1/b —1]a 


Means of Two Numbers 
Arithmetic mean = A = }(a + b) 
Geometric mean = G = 4/(ab) 
Harmonic mean = H = 2ab/(a + b) = G?/A 


Sums of Important Series 
Sum of the first n natural numbers = $n(n + 1) 
Sum of the squares of the first n natural numbers 


= 4n(n + 1)(2n + 1) 
Sum of the cubes of the first natural numbers = [4n(n + 1)]? 


1.9. COMBINATORIAL FORMULAE 
Any arrangement of n objects in a definite order is called a 
permutation. 
The number of permutations of n objects amongst themselves is 
n(n —1)n—2)...21=n! 
The number of permutations of n objects r at a time is 
"P, =n!/(n —r)! 


The number of combinations of n objects taken r at a time is the 
number of selections, each of r objects, without reference to 
their order within a selection, and this number is denoted by 


SC, or (Ty 
ar 


(") = n! a, MOT 2). .. nF +2) 
rr) r\n—P! r! 


The number of ways in which n objects may be divided into k 
groups containing respectively m, M2, Ms,. . .,M, Objects such 
that m +ng+...+m =nis 


n!/(ny! Mal oc. ny!) 
1.10. BINOMIAL AND MULTINOMIAL THEOREMS 


Binomial Theorem 


ctor=|((tpe (sere (Soret. 


+(")eart...+ ("Jarl 


(n is a positive integer.) (See also Binomial Series, p. 44.) 


The values of the binomial coefficients for different values of n 
are: 


Value of n Binomial coefficients 
1 1 1 
2 1 2 1 
3 1 3 3 1 
4 1 4 6 4 1 
5 1 5 0: 100 1 


etc. etc. 


Note that the binomial coefficients for any value of n may be 
obtained by summing consecutive pairs of coefficients in the 
preceding line. 


Multinomial Theorem 
(X1 + Xa +... + x)" 


n! 
= > xy". tee 
my! Mal. . . my! 


ke 
where n = 2 n, and each n, takes integral values from 0 to n. 
i=1 


1.11. DETERMINANTS 


Second order determinant | a; 5, | has the value (a,b2 — ab) 
a2 be 
Third order determinant | a, 6; c, | has the value 
a2 be Ce aA, + b,B, + eC, 
a3 bs C3 
where Ay = ‘be Ca B, =] Ce Ae Ci =] a2 be 
bs C3 C3 As a3 bs 


The solutions of the simultaneous equations a,x + biy = cy 


a2x + bey = C2 
are given by 
x 1 
Cy by a Cy a b, 
C2 be a2 C2 ag be 


and the solutions of a,x + biy + cz = d, are given by 
2X + bay + Coz = dz 
asx + bsy + C3z = dg 


b a y a Zz Re? Pe We 
d, by cy a a cy a bh da a b; C1z 
dz be C2 a2 dz C2 az bo dae az bz C26 
ds bs C3 a3 ds C3 az bs ds a3 bs C3! 
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MENSURATION AND 
TRIGONOMETRY 


2.1. ANGLE 


One radian is the angle subtended at the centre of a circle by an 
arc equal in length to the radius r. 


27 radians = 360° 

7/180 radians = 1° 

1 radian = 180°/7 = 57:3° 

Length of arc = s = ré (@ in radians) 


2.2. TRIGONOMETRICAL RATIOS 


Let 0 be the angle through which 
the radius OA would turn to 
reach the position OP, ie. 6 = 
ZPON. 

(Convention: counter-clockwise 
rotation gives a positive angle 
and clockwise rotation gives a 
negative angle.) 


ZOPN = 90° 


ON and PN are directed lines in 
accordance with the convention 


for Cartesian coordinates so that 
ON is positive in the Ist and 4th quadrants 
and negative in the 2nd and 3rd 
NP is positive in the 1st and 2nd quadrants 
and negative in the 3rd and 4th 


DEFINITIONS 
sin 0 = NP/OP cos 6 = ON/OP tan 6 = NP/ON 


Algebraic signs of the above ratios in the four quadrants 
+ + _ + _ + 


sin cos tan 


Further Definitions and Identities 


cosec @ = 1/sin 0 sec 0 = I/cos 0 cot 6 = 1/tan 0 
sin~*x is the angle 0(— 42 < @ < 47) such that sin 90 = x 
cos~*x is the angle 6(0 < @ <7) such that cos 6 = x 
tan~+x is the angle 0(— $7 < @ < 47) such that tan 6 = x 


Identities 
sin? 0 + cos? 0 = | 1 + tan? 0 = sec? 60 
1 + cot? 6 = cosec? 6 tan 0 = sin 6/cos 6 


Complementary Angles A + B = 90° 
sin A=cosB cosA=sinB tan A=cotB 


Supplementary Angles A + B= 180° 
sin A = sin B cos A = —cos B tan A = —tan B 


2.3. MENSURATION 
Plane Areas 


Parallelogram Sides a, 6 and included angle «. Perpendicular 
distance between parallel sides of length a = h = b sin « 


Area = ah = ab sin « 


Trapezium Area = 4 (Sum of parallel sides) x Perpendicular 
distance between them 


Circle Radius r, diameter d = 2r. 
Circumference = 27r = ad = Area = mr? = }nd? 
If an arc subtends an angle @ radians at the centre 


Area of sector = $r?0 
Area of segment = $r? (0 — sin 6) 


Ellipse Semi-axes a, b; area = mab. 


Solids 


Notation: V = volume, C = curved surface area, S = total 
surface area. 


Solid Circular Cylinder 
h = height, R = radius 
VeanR*h C=27Rh S = 27Rh + 27R? 


Hollow Circular Cylinder 
r = inner radius 
V = (R? — rh 


Pyramid and Cone 
V = 3 (Area of base) x Perpendicular height 


Right Circular Cone 
h = height, R = radius of base, / = slant height = /(R? + h?) 
V = 4rR*h C= wRi 
S = 7R? + 7RI/ 
Volume of frustum = 47h (R? + Rr + r?) 
where R, r are radii of the ends and h is the height of the frustum. 
Sphere Radius R V = $7R° C = 4nR? 
Volume of a zone = @7h [3(R? + r?) + h?] 


where R, r are radii of the ends and A is the height of the zone. 


2.4. SOLUTION OF TRIANGLES 


Triangle ABC: sides a, b, c. Perpendicular height above the 
base a is h. Angle A is opposite side a, etc. R, r are the radii of 
the circumscribed and inscribed circles of triangle ABC. 


Semi-perimeter = s = 3a + 6 +c) 
Area = $ah = fabsinC = \/[s(s — a)(s — b)(s — od) ] = rs 
If C = 90°, c? = a? + b? (Pythagoras’ Theorem) 

Sine formula: a/sin A = b/sin B = c/sin C = 2R 


Cosine formulae: a? = b? + c? — 2be cos A ete. 
a=bcosC +ccosB ete. 


sin $A = fi BS = 2) ie ae =f cos 4A = Y pil = ) 


(s — b)(s —c) 


tan 3A = Gee 


Area of a cyclic quadrilateral (sides a, b, c, d) 
VI(s—a)(s—b)(s—c)(s—d)] wheres =}}(a+b+c+d) 
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2.5. COMPOUND ANGLES 


sin (A + B) = sin A cos B + cos A sin B 
cos (A + B) = cos A cos B F sin A sin B 

tan A + tan B 
tenis eae Fe ign A tan B 
sin 2A = 2 sin A cos A 
cos 2A = cos?A — sin?A = 2 cos?A — 1 = 1 —2 sin?A 
sin 3A = 3 sin A — 4 sin® A 
cos 3A = 4cos*A — 3 cos A 
sin C + sin D = 2 sin 3(C + D) cos 4(C — D) 
sin C — sin D = 2 cos 4(C + D) sin #(C — D) 
cos C + cos D = 2 cos (C + D) cos #(C — D) 
cos C — cos D = —2 sin 3(C + D) sin 4(C — D) 


2.6. SPHERICAL TRIGONOMETRY 


Spherical Triangles 

The section of a sphere by a plane 
through its centre is called a great 
circle. A spherical triangle is the figure 
bounded by three great circles. 

The sides BC, CA, AB of the 
spherical triangle ABC are measured 
by the angles a, b, c, which they subtend 
at the centre O, The angles A, B, C are 
the angles between the planes of the 
great circles. 


Spherical excess = E=A+B+C—a 
Area of spherical triangle ABC = (A + B + C — n)r? = Er? 
Sine formula: sin A/sin a = sin B/sin b = sin C/sin c 
Cosine formulae 
cosa =cosbcosc + sinbsinccosA ete. 


Cotangent formulae 
cot b sinc = cos ccos A +sin A cot B etc. 


cos A + cos Bcos C 


cosa = a ain ae etc. 
tan (A + B) ee iC 
tan (4 — B) = see cot 3C 
tan 4a + 6) = eR tan 4c 
iia) — 2214 - |e 


sin 4(A + B) 
Ifa+b+c=2s 


she /—-=ee-e (s — b) sin (s — c) 


sin 6 sin c 


sin s sin (s — a) 
sin b sine 


ahd = Ri eet 2 den (s — b) sin(s — c) 


sin s sin(s — a) 
IfA+B+C=2S 
— cos Scos (S — A) 


cos $A = 


eae sin B sin C 
See cos (S — B) cos (S — C) 


sin B sin C 

— cos Scos (S — A) 
t SS —— eee 
an ae lane eS C) 
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Delambre’s Analogies 
sin 3(A + B) cos 4c = cos (a — b) cos $C 
sin 4(A — B) sin $c = sin (a — b) cos $C 
cos 3(A + B) cos 4c = cos #(a + b) sin $C 
cos 3(A — B) sin $c = sin 4(a + 5) sin $C 


Section 3 


COORDINATE 
GEOMETRY 


3.1. DIVISION OF LINE IN GIVEN RATIO 

P (x, y) divides the straight line AB in ratio m:n, i.e. 
AP/PB = m/n 

With A (x, y1) and B (xa, ya), 


x = Mats = 1+ mya 
m+n m+n 


[For external division, the sign of m becomes negative.] 


3.2. EQUATIONS OF A STRAIGHT LINE 


Standard form m = tan 6 = gradient or slope, c = intercept 
on y axis 

y=m-+e 
Perpendicular form p =\ength of perpendicular from the 
origin O to the line and « is the inclination of the perpendicular 
to Ox 


xcosa-+ysinan=p 
Intercept form aa, b are the intercepts on Ox, Oy respectively 


< pe 
ze ; 
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General form 
Ax + By+C=0 
y—ji =m(x—-x,) [Line of slope m passing through 
(x1, yi] 
Bi Sit A 


yr-n= (222) (x — x1) [Line passing through (x1, y1) 
we M1 and (x2, y2)] 


Angle between two lines y = mx + c, and y = mx + Co is 


m, —m 
NE ea 
1 + m, Me 
Two lines are perpendicular if mymz = —1. 


Laws from Plotted Results. Form of law modified to give 
straight line graph. 


Law Other Quantity Quantity Slope Inter- 
suspected formof to be to be gives: cept 
law plotted _ plotted gives: 
along along 
horizon- vertical 
tal 
y=ax*+b —__— x? y a b 
xy = ax 
+ by yobylxta ylx y b a 
y = ax" logioy = logiox logioy un logioa 
nlog x 
+logioca 
y = ace logioy = x logioy 0:4343k logioa 
0:4343kx 
+ logioa 
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3.3. CIRCLE 
(i) x? + y? =a? Circle centre origin, radius a 
(ii) (x — h)? + (y — k)? =a? = Circle centre (A, k), radius a 


(iii) x? + y? + 2ex +2fy +ce¢=0 General equation of 
J & y q 
circle 


Tangent at (x1 y:) on circle (i): 

xX1 + yy, = a? 

Tangent at (x1, yi) on circle (ii): 

(x — Arr — A) + (y — kK) (1. — k) = ? 
Tangent at (x1, y1) on circle (iii): 

M+ teetxua)t+fyty)+e=0 


3.4. CONIC SECTIONS 


A conic section is the locus of a point which moves in a plane 
such that the ratio e (called the eccentricity) of its distance from 
a fixed point (the focus) to its distance from a fixed straight line 
(the directrix) is constant. 

The locus is an ellipse, a parabola or a hyperbola according 
as 


e<=>1 


Polar equation of a conie referred to its focus where / = semi- 
latus rectum 


: =1+ecos 0 
, 
(a) Ellipse 
Ellipse, centre origin, semi-axes, a, b, e < 1 
x? 
b? = qa? (1 — e?) l= b?/a 
16 


Tangent to above ellipse at (x;, y:) on it 


xX 
athe 


(6) Parabola 
y* =4ax Parabola, origin at vertex, Ox axis of symmetry 
yyi = 2a(x + x1) Tangent at (x1, y1) on the parabola 


y=mx + — Tangent of slope m for any value of m 


Parabolic arc between two points at the same level. Span /, 
height s at the mid-point referred to origin at one end: 


y= xl — x) 


Area of segment of parabola = $ s/ 


(c) Hyperbola 

x? y? 

ae = Hyperbola, centre origin, semi-axes a, b 
pal ee 


7m cia. Tangent at (x;, yi) on the hyperbola 


~ + % = 0 Equations of the asymptotes 


b? = a?(e7—1) (e>1) 


xy =a’ Rectangular hyperbola, asymptotes are the coordinate 
axes 


yiX + X1y = 2a? Tangent at (x1, y1) on rectangular hyperbola 


(d) Focal Distances 


The ellipse and the hyperbola have two foci and two directrices 
whilst the parabola has only one focus and one directrix. A 
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focal distance is the distance of a point on the conic from a 
focus. 


Ellipse: the two focal distances are a + ex. Sum = 2a. 
Hyperbola: the two focal distances are ex + a. Difference = 2a. 


3.5. GEOMETRY OF SOLIDS 


(a) Point 


Let (x, y, z) and (r, 0, 4) be the Cartesian and Polar coordinates 
of a point P, 


x=rsinOcos¢ y=rsin@sin¢ z=rcos 0 


where r = OP = \/(x? + y? + 2?) 
6 =co-latitude and cos 0 = z/4/(x? + y? + 2?) 
¢ = longitude and tan ¢ = y/x 


(6) Straight Line 
Distance between two points P (x1, Yay 21) and Q (x2, ya, Z2) 
d = V[(%1 — X2)? + (1 — yo)? + (21 — 22)] 


The Direction Cosines of a straight line PQ are the cosines of 
the angles between PQ and Ox, Oy, Oz and are denoted res- 
pectively by /, m, n. 


Xo = X — Za Z 
[et Gh ag: oi pp ee eee 8 ae 22 2 


d d d 
P+m+n=1 


The angle @ between two straight lines whose direction cosines 
are (/,, m1, 7) and (lz, me, nz) is given by 


cos 6 = hl, + mym2 + nyNn2 
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Equation of the straight line passing through the point (a, b, c) 
having direction cosines (/, m, n) is 


—_— SO 


(x, y, z) = 0 represents the equation of a surface 
(c) Plane—linear in x, y, z 
General equation of a plane 
ax + by +cz+d=0 
Intercept form: plane cuts the axes at distances a, b, c from O 
ero ce 


Perpendicular form: perpendicular from O to the plane has a 
length / and direction cosines (/, m, 7) 


Ix + my + nz =p 
Length of perpendicular from (x1, yi, 21) to the above plane 
Ixy + my; + nz, — p 


(d) Surfaces of the Second Degree 


Surfaces of revolution 


x? + y? =a? Right circular cylinder, radius a 
x? + y? = a®z? Right circular cone, vertex origin, axis Oz 
x? + y? +z? =a? Sphere, centre O, radius a 


ee ae : 
a “f+ pa + rt 1 Ellipsoid, semi-axes a, b, c 
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Section 4 
NON-DECIMAL ARITHMETIC 


4.1. ARITHMETIC BASES 


Decimal Notation: numbers are expressed as multiples of powers 
of 10. 


10 is called the base. The symbols used are 1, 2,. . ., 9. 
293 means 2 x 107 +9 x 10 +3 


Any other positive integer may be used as the base of a number 
system, e.g. in octal arithmetic, 537 means 5 x 87 + 3 x 8 +7. 
Decimal equivalent: 351 

More generally, in arithmetic base r, the number represented by 


Gn An-1.- . «42QA, a 
where each of the digits ao, a1, d2,. « +» @n-1)@n < (an # 0), is 
Ant” + Anaat™ 2 +... + ar? + aur + a 


Binary arithmetic (base 2), using the two digits 0, 1 only, has 
become very important because it is widely used in digital 
computers. 


4.2. CHANGE OF BASE 


N is an integral number expressed in any base. It is required to 
express N in base r. 
Let N = agr® + Qn-ar™ + +... tar tao. 
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Divide N by r so that N/r = agr*~* + Qn-ar"-2 +... +41: 
remainder do. 
Divide the quotient by r. The next remainder is a, and so on. 
The successive remainders are the required digits ao, a1, . . ., Qn. 
For example, by successive divisions by 2, decimal 247 becomes 
11110111 in binary. Binary numbers are most easily converted 
to decimal by reversing the process 
iwi «68M 1 1 
Be ts 30 61 123 247 


Multiply the left-hand digit by 2, add the next binary digit, and 
continue this process. 


4.3. FRACTIONS 


In decimal arithmetic, 9-367 means 9 + = ig cape oh 


08 * 10" 
Similarly in binary arithmetic, 1.101 means 1 + 5+ +5. 
In arithmetic base r, the fraction 0.5; b2b3. . . where 5, ba, b3, 

. . are integers, each less than r, of which one or more may be 
zero, is equivalent to 


ee. be 
S+at+Gt... 


To express any fraction F in the scale of r, bi, b2, bs, . . . must 
be found such that 
by, be bs 
F=—+—-+—+... 
: 7 a ~ 3 - 


Multiplying byr = Fr=6,+ “2 +23 Wis 06 


so that b, is the integral part of Fr, the fractional part F, being 
given by 


and bz is the integral part of Fir and so on. 
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For example, to express 11/16 (decimal) as a binary fraction: 
X2=1+8 FxX2=042 Ex2=144 
4$x2=1+0 


11/16 (decimal) = 0.1011 (binary) 
4.4. THEOREM 


In the scale of r, the sum of the digits of any integer when 
divided by (r — 1) gives the same remainder as when the 
integer itself is divided by (r — 1). 


Corollary: in the decimal scale (r = 10), a number is divisible 
by 9 if the sum of its digits is divisible by 9. 
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Section 5 
COMPLEX NUMBERS 


5.1. COMPLEX NUMBERS 


If a and b are real numbers and i = +«/ —1 (i.e. i? = —1), ib is 
known as an imaginary number and z = a + ib is called a 
complex number and Z = a — ib is called the conjugate of z. 

If zi =a+ib and z2 =c + id, their sum, difference and 
product are 


4+22=(@aic)+ib+d). 
232g = (ac — bd) + i(ad + bc) 


zZ = q? + 5? 

Z,_a+ib_(a+ ib(c—id) y 

Z2 c+id (c+id)(c — id) P(x, y) 

_ (ac + bd) + i(bc — ad) . (r, @) 
c? + d? 

If z,; = zo, thena=c,b=d. 8 x 


O 


5.2. ARGAND DIAGRAM Argand diagram 


There are two geometrical representations of the complex 
number z = x + iy: 

(a) the point _P(x, y) (6) the vector OP 
| 2=x+iy=r(cos 6 + isin 6) 
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DEFINITIONS 

(a) Absolute Value or Modulus of z = |z| = |x + iy| 
= length of OP =r = (x? + y?)1/2 

(6) Argument or Amplitude of z = arg z or amp z = 0 
= tan~? (y/x) where —7 < 0 <a. 


, 7 | | I 
ol tea gee 
Se a (ae 
=(1 ata 2) $716 a ag 


= cos @ + isin 0 
Therefore e-'® =cos@—isin@ so that 
cos 6 = $e? + e-'%) = cosh 10 
isin 6 = }(e — e~%) = sinh i6 
Therefore i tan 6 = tanh i0 
Also z=re'® and Z=re~‘ 
2122 = r,(cos 6; + isin 0,). re(cos 02 + i sin 62) 
= r,e'91 , ree! 
= ryree'(1 + 92) = ryro(cos 0, + O02 + isin 6, + 02) 
arg 2172 = 0; + 02 = arg 2 + arg z2 
[2122] = rire = |Z;| . [zal 


5.3. DE MOIVRE’S THEOREM 

{r(cos 0 + isin 6)}" = r"(cos nO + isin nO) for all values of n 
nth roots of z 

z=x+iy =r(cos 0 + isin 6) 


= r{cos (0 + 2kz) + isin (0 + 2km)} = retO+2km) 
where kK = 0, 1,2,.. ., (2 — 1). 


B/ z= "/(x +iy) = ri {cos (Cte) + isin (- adil} 
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nth roots of unity 
1 = cos 0 + isin0 therefore 


*/1 = cos (=) + isin (==) er aay ee 


‘-. 
Cube roots of unity (” = 3) 


24/1 = cos (=) + isin (=) (k =0, 1,2) 
The three roots are 1, —4$ + i\/3/2. 
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Section 6 
MODERN ALGEBRA 


6.1. SET 


A set is a well-defined collection of objects. ‘‘Well-defined” 
means that it may readily be decided whether an object is a 
member of the set or not, e.g. 


(i) The set of all counties in the U.K. The elements of this set 
may be listed: {Aberdeenshire, Anglesey, .. . Yorkshire}. 
(ii) The set of all positive even integers less than 16: {2, 4, 6, 


* "9 


A set may contain an infinite number of elements. 


6.2. NOTATION 
x € A means x is an element of the set A. 
(x ¢ A is the negative statement.) 


A = {X1, X2,. . ., Xn} Means x1, X2,. . ., X, are all the elements 
of set A. 


A = {x, P} means 4 is set of elements x satisfying a property P. 
3 means “‘there exist(s)”’ 
For two statements p and q, 
P > q means “‘p implies q” or “if p then q”’. 
p <q means p and q are equivalent statements or “p if and 
only if g”’. 


26 


A = B means two sets A and B are equal, i.e. A and B contain 
precisely the same elements or x € A <> x € B. 

A ~ B means sets A and B are equivalent, i.e. there is one- 
to-one correspondence between their elements. 

A © Bmeans A is a sub-set of B, ie. every element of A is 
also an element of B. 
AC B<+(xEeA> xEB) 
A=B<Ac BandBCA 


6.3. OPERATIONS ON SETS 


A U B denotes the union of two sets A and B. It is the set of 
those elements which belong to either A or B or both 
A and B. 
ie. xE€AorxEBoxEAUB 

Aq B denotes the intersection of two sets A and B. It is the set 
of those elements which are members of both A and B. 
ie. xE€ AandxEBoxEANB 

The universal set V is the set of all possible objects in a given 

problem. The empty or null set ¢ is the set having no elements. 

If two sets A and B have no elements incommon, AN B = ¢. 


6.4. VENN DIAGRAMS 


Let the universal set V be represented s' 

by a rectangle and let any subset S Vv 
of V be represented by.a circle 
(shaded) drawn inside the rectangle. 
The set S’ (unshaded) of all elements 
not in Sis called the complement of 
Sand Su S’ = V. 


Aja Bis represented by the doubly 
cross-hatched area. 


A VU B is represented by the whole 
shaded area 


y 3) 


The following important results can easily be verified by Venn 

diagrams: 

i) (ANB)ANC=AN(BNQG) age: 
(AUB)UC=AU(BU oI Associative Laws 

fii) AU(BNQ=(AVBAN(AVUC) 
AN(BUC)=(ANB)U(A NC) 


(iii) 4 pie 4 i os 7 7) De Morgan’s Laws 


Distributive Laws 
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Section 7 
MATRICES 


7.1. THE ORDER OF A MATRIX 


A matrix is an ordered array of numbers or other elements. 
A matrix having m rows and n columns is said to be of order 
m Xn. 


411 Ae As 2\ is a column matrix; 
421 A222 Gaz 1} order 3 x 1. 
431 32 433 3 


G41 Gag Gag 
is a4 X 3 matrix. 
(—2 1 O 7)isarow matrix; order 1 x 4. 


A matrix having ” rows and n columns is called a square matrix 
of order n. Matrices will be denoted by heavy capital letters 
a 


7.2. THE MULTIPLICATION OF MATRICES 


Suppose A is a 3 X 5 matrix and Ba 5 x 4 matrix. Then, for 
example, the element cos of the product A x B is obtained by 
the multiplication of the elements of the second row of A by the 
corresponding elements of the third column of B, adding the 
products so that 


C23 = Aaibisg + deabe3 + Aagb33 + Aeabaa + Aosbss 


5 
= 2 Aadys 
k=1 
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In general 


5 
Cy = De aucdis 
k=1 


Compatibility 
The above product rule is applicable only if the number of 
columns of A equals the number of rows of B. Matrices which 
satisfy this criterion are said to be compatible for multiplication ; 
otherwise, their product is not defined. 

In general, a matrix A (m x n) is compatible with a matrix 
B (n X p) and A x Bis a matrix C (m xX p). 

A matrix A (m x n) is compatible with a matrix B (n x m) 
for multiplication in either order but A. B # B. A in general. 


7.3. UNIT MATRIX: NULL OR ZERO MATRIX 
(se - = 969 
S Ayre d c d ec dj\O0O 1 
If any (2 X 2) matrix is pre- or post-multiplied by I = ( 0 4 


it is unchanged. 
I is called the unit matrix of order 2. 


100 
(0 1 0) is a unit matrix of order 3, and so on to higher 
00 1 


orders. 
A zero matrix has all its elements zero. 


7.4. EQUALITY; ADDITION; MULTIPLICATION BY 
A SCALAR 


(a) Two matrices are equal if and only if they are of the same 
order and their corresponding elements are equal. 
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(6) The sum of two matrices of the same order is defined to be 
a matrix whose elements are the sums of the corresponding 
elements of these two matrices. 

(c) The product of a matrix by a number (scalar) k is defined to 
be the matrix derived from the first by multiplying each of its 
elements by K, e.g. 


te 6 _ [ka kb 
a= (2 ) then kA = (12 ta) 


7.5. DETERMINANT OF A SQUARE MATRIX 


Matrix A (n x n) Determinant of matrix A = |A| 
M1 Ayo...Ain Q@y1 Ayg. . - Ain 
G21 22 aon G21 Aeo2 Qan 
Ani Ane ann Qn1 an2 Ann 


A rule assigns a value to the determinant [A] of order n. 


Let a, be the square matrix of order (n — 1) obtained by 
deleting the ith row and the jth column of A. The determinant 
|et4s| is called the minor of ay in [A]. The cofactor Ay of ay is 
defined by 


Ay = (—1)'*! x value of [e;5| 
The value assigned to the determinant |A| is 
n 
|A] = & ay Ay 
j=t 
(This formula will give the same value of |A| irrespective of the 
particular row in terms of which it is evaluated.) 


If the products obtained by multiplying the elements of any 
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row by the corresponding cofactors of another row are added, 
the result is zero. 


n 
py ay, Ay, =O Ci #j) 
k=1 


oe bss 7) 
(2x2 Q@21 422 
a a 
JA] =] O22 22 1 = ay Air + 12412 = 411422 + A12 (—aa1) 
421 422 


441 Qa Aig 
A=1|4@21 Qe2 dag 
(3x3) \ag, ds2 ss 

41 G2 aig 
421 422 423 
431 432 433 


= 411A11 + Q12A12 + @13A13 


JA] = 


G21 422 
431 Ase 


Q21 Qe3 
Q31 433 


G22 Gaz 
= a4 — Qi2 + ais 
Q32 §=az3 


A matrix whose determinant is zero is a singular matrix. 


One of the rules for multiplying two determinants of the same 
order is precisely the rule (Sec. 7.2) for multiplying two square 
matrices of the same order. Hence if A, B are square matrices 
of the same order 


|AB| = |A]. |B] and |AB| = |BA| 
even though, in general, AB 4 BA. 
7.6. INVERSE MATRIX 


The inverse of a square matrix A is A~+ of the same order such 
that 


A“?*A = AA“? = 
32 


xe EE b t= 7 ( d x 
IA = (2 4), then A “ail a 


A singular matrix has no inverse. 


7.7. MATRIX SOLUTION OF LINEAR EQUATIONS 
Equations ax+by=p and cx+dy=q 


: aa B\ix\ ip 
Matrix form: : 4 (>) = (F) 
Solution: (*) pga (pa sity ") 
A |A] \99 — ¢p 
7.8. THE TRANSPOSE OF A MATRIX 


If the rows and columns of matrix A (m xX n) are interchanged, 
the resulting matrix A’ ( x m) is said to be the transpose of A. 
A square matrix A is symmetric if A = A’ and skew symmetric if 
A’ = —A. 

If A = A’ (symmetric), the elements a,; of A are such that 

ayy = Ay. 

If A’ = —A (skew-symmetric) a,; = —ay,, (i #/), whilst ay, = 0. 


7.9. ADJUGATE OR ADJOINT MATRIX 


If A is a square matrix, the adjugate matrix adj A is formed by 
replacing each element of A by its cofactor in |A| and transposing, 
so that 

adj A= (Ay;)’ = (A;;) 
where (A,,) denotes a matrix of order n having elements Ajj. 


4341 G19...Ain Ass Am. ..An 
Q21 A422 dan Ayo Ane Ano 
= then 
oe : adj A = 
QAn1 Ane ann Ain Aon Aan 
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7.10. INVERSE OF A SQUARE MATRIX 


If A is a non-singular matrix of order n, 


1 4 
A-i = —adjA 
[A] 
7.11. SOLUTION OF LINEAR EQUATIONS 
(@) yX1 + ayaXo t+... + AnxXn =O, G@=1,2,.. ., n) 
are n simultaneous linear equations in n unknowns x. 
Matrix form AX = B 


where A is a square matrix of the coefficients a,; and X, B are 
column matrices of the x, and b,. 


If A is non-singular, X = A~*B where A~* is given in 7.10. 


(6) Method of Triangular Decomposition illustrated for three 
unknowns. 


Matrix form AX =B (Ais the 3 x 3 matrix of a,,) 


1 6.8 W141 U2 Us 
Find V = (x 1 0) and U = (0 U22 in) 


V31 UVsq 1 0 0 Usa 
such that VU = A. 
From the equations, Uy, = 411, Ug = QAi2, - - +» Vai Ui1 = 421, 
a. oa ase 


Let VY = B(Y is a3 X 1 matrix) whence y,; = 5; and 
Vai yi + Yo = be Vai Vi + Vs2Y2 + Ys = bs 


from which Y is found. 

AX=B VUX=VY UX=Y 
whence 3%1 + Ui2X2 + UisXs = yi and 
Uz2X2 + UasXs = Yo = UaaX3 = a 

from which X is found. 


34 


Section 8 
VECTORS 


8.1. COLUMN VECTOR 
—> 
Position vector of P(a, b) = OP = 1. 


r, may be represented by a column matrix or column vector (7) : 
Addition of Two Vectors. Let r2 = (4) , then 


mtr = (5) + ( ‘) = ( ri H 4) (Parallelogram Law) 


y R(atc, b+d) 


Base Vectors 


i= ( 0) and j = (9) represent unit vectors along Ox, Oy 


respectively. 
r= (3) =a(}) +6(%) = ai + bj 
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8.2. LINEAR TRANSFORMATIONS IN A PLANE 


A column vector & ) is transformed into another column vector 
(>) on pre-multiplying the former by a 2 xX 2 matrix A 


(transformation matrix). 
ei wk te) ay (* 
(;) hd (3 4) 7 


ie. r’ = Ar and x’ = ax + by, y = cx + dy 


The base vectors i, j transform into i’, j/ such that 


a= (2 4) (8)=(2) = as a= (8) =1 


i.e. i’, j/ are the columns of A. Thus if the images of the points 
(i, 0) and (0, 1), under a transformation, are known, the 
transformation is known. 


8.3. AFFINE TRANSFORMATIONS 


ee Fae + by) a ee s 
(7) = (Sth) = * (0) + (4) = 491 
The point (x, 0) transforms into (ax, cx), ie. x (a, c). 
The point (0, y) transforms into (by, dy), i.e. y(b, d) 


A grid of unit squares in the (x, y) plane will transform into a 
grid of parallelograms. The transformation is such that parallel 
lines remain parallel and the ratios of lengths in the same 
direction are unchanged. It is called an affine transformation, 

Exceptionally, if the matrix A is singular, i.e. if |A| = 0, then 
alc = b/d or a = kb, c = kd so that x’ = (b/d)y’ and all points 
in the (x, y) plane are transformed into points on one straight 
line. 
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8.4. BASIC TRANSFORMATIONS 


Type Matrix Algebraic 
Equivalent 
Reflection in Ox (4 4 x hg 
0 —1 yory 
F —-1 0 x’ = —x 
Reflection on Oy ( 0 i) yey 
Se ee ost 0 1 x =y 
Reflection in line y = x (i a) y=x 
erin 0 -1) x =-y_ 
Reflection in line y x * 1 0) y= —% 
woe 0 —!1 x= -y 
Quarter turn about origin ( 1 4) y =x 
me —1 0\ x’ =-x 
Half turn about origin ( 0 wae yo-y 
¥ «os a | x =y 
Three-quarters turn about origin ( an i) y wee 
r & x’ =kx 
Enlargement ( 0 i) y =k 
, .. oe xX =x+k 
Shearing parallel to Ox ( 0 ) eae y 
Rotation through angle « (= oan | 
sina cosa 


x’ =x cosa — ysin « 
y =xsin a + ycos « 


Reflection in line y = x tan « Res 2% sin ao 


sin 2a —cos 2a 
4 


x’ = x cos 2a + y sin 2a 
y’ = x sin 2x — y cos 2a 
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8.5. EIGENVECTORS 


Consider the transformation 


GAG)“ G) 
"y y e d/\y 
where the transformation matrix A is non-singular. Vectors 


(3 whose directions are unchanged by the transformation are 


called eigenvectors, i.e. ( a = ( kx ) 
" y ky 


We require solutions of (¢ 4) (*) BS k(*) 
c d}\y y 


ie. of (A — aD(*) =) 


The eigenvalues of A are the values of k obtained from the 
determinantal equation 


|A — kl|=0 
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Section 9 
DIFFERENTIAL CALCULUS 


9.1. LIMITS 


A function f(x) is said to approach the limit ZL as x approaches c, 
ie. lim f(x) = L if, given any positive number «, there is a 
positive number 6 such that 

lf) —L|<e when0 <|x —c| <d 
Some important limits 
sin 0 _ im t20 0 


9.2, DIFFERENTIATION 
Let y = f(x) and let dy, dx be corresponding increments in y,x 
then the differential coefficient or derivative of y with respect to x 
is given by 
Oe ae tim 2 + oo) SIO) 
dx 7-0 OX 82-0 ox 
= gradient of the curve y = f(x) at the point (x, y) on it. 
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Other symbols for this differentia! coefficient are y’, Dy, and 


f'). 


Rules for Differentiation 


Let u, v, w, etc. be functions of x. 


y = cu (c = const) 


yroutov+wet. 


y =u 
=H 
Ie 
y = uw 
M Soe 
y=f&) 


where u = ¢$(x) 


If y and x are both functions of f, 


If x,y are distances and ¢ is time, 


x,y) are velocities. 
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dy _ du 
dx dx 
dy du, dv , dw 
ay ere ete sh. - 
dy_ adv du 
de ae ie 
o=5 (oS du _ u®) 
dx dx Oe 
2 yoy (14H 4,1 do 4 1 de) 
dx udx vdx wdx 
AY yn dt 
ae ae 
dy __dy du_df du 
dx du dx Me dx 
dy _ dyldt 
. dx|/dt 


&. S ; (often denoted by 


Table of Important Differential Coefficients 


y dyldx 
x” nx” - 1 
eks kekz 
sin (ax + b) a cos (ax + b) 
cos (ax + b) — asin (ax + 5) 
log. (ax + 5) al(ax + b) 
tan x sec? x 
cot x — cosec? x 
sin=+ x/a 1/\/(a? — x?) 
cos~* x/a — 1//(a? — x?) 
tan~1 x/a al(x? + a?) 
sinh x = 4(e* —e-*) cosh x = ke? + e-*) 
cosh x sinh x 
e* pears e-* 
tanh x sech? x = 1/cosh* x 


e* -++ e7-? 
9.3. NEWTON’S METHOD OF APPROXIMATION 


If xo is an approximation to a root of f(x) =0, a better 
approximation x, is given by 


eA f(x) 


XO = X11 
so that Xi = Xo — f(xo)/f’(xo). 
9.4. SUCCESSIVE DIFFERENTIATION 


If dy/dx is differentiated with respect to x, the second deriva- 
tive d?y/dx? is obtained. Successive differentiations produce 
derivatives of higher orders. 
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Notation 


ty = fo, 2 (2) -2=4yeo=s'@ 
(ry 7 ee 
uel ao 


nth derivative = 2 = Dry = f(x) 
Leibniz’ Formula for the nth derivative of a product. 


D"(uv) = Du. o+(7 ) p»- “a. Do + (3 )p- 2u, D?v+. 
n * n . 
+(7)> u.Dot+...+(,"4)Du. dD" ty tu. Dv 


9.5. MAXIMUM AND MINIMUM VALUES 
If oa = 0 at x = a, then 


= f(x) has a minimum at x = a if d*y/dx? > 0 at x =a 
chek s mesinom at x =o i dy /dx?< Oatx=a 


If ‘Pde = 0 at x = a, the graph of y = f(x) hasa point of 
inflection at x = a. 


9.6. CURVATURE 


y = f(x) graph APB 
PT, PNare tangent and normal 
at P respectively. 


dy/dx = tany 
= slope of tangent 
ae ae 
Subtangent QR dyldx 
a yiy 
Subnormal QN = =-»2 = yy’ 
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Length of tangent = PR = y(1 + y%)*/?/y’ 

Length of normal = PN = y(1 + y”%)*? 

Let s be the length of the arc AP measured from a fixed point A. 
dy 


Curvature = «x = ne yA + 2-82 
Radius of curvature = 1/« 


If the curve is nearly horizontal at all points as in the case of 
horizontal beams, dy/dx is small compared with unity and the 
curvature is d?y/dx* approximately. 


9.7. PARTIAL DIFFERENTIATION 


Itz = f(x,y), de = 2 de + = = dy = Lact Za 


dx, dy, dz are called sm itale 
dz_ ozdx , dzdy 


If x, y, z are all functions of f, = ee alee § 


If u = f (x1, X2, Xs, . + +, Xn), then 
du = dey + dita + diy + +e dts 


Implicit ape If y is an implicit function of x as given by 
f(, y) = 
Cini ¥en 6 .2=- Ze 
ox oy 


9.8. INFINITE SERIES 


A series may contain either a finite or an infinite number of 
terms. Let S, denote the sum of terms of a series. The series 
is said to be convergent or divergent according as the lim S, is 
a finite quantity or is + ©. -jiaiied 
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Binomial Series With reference to 1.10, if n is a positive integer, 


atone (8) +(—)e# Ghee ()e 
5 . +( Ghee 
= 1+ nx +5 mn — x? +. ..t+ 


< n(n — 1). ~. @—r+tix t+... +x" 


which is’a finite series. 
If n is not a positive integer 


(La = 1 4 nx + oan — Dat +... + 


ann — 1). art xt +... ad inf, 


provided |x| <1, ie. the infinite series converges to (1 + x)" 
provided |x| < 1, otherwise it is divergent. 
Approximations. If x, y, z, u, v are small compared with 1, 
(1 +x)"=1+nx foraliln 
e.g. /(1 + x) = 1 + 3x 


(1 + xl + yl + 2) 
(1 + wu) + v) 


sS1l+x+ytz—-u-—v 

Exponential Series 

eh mel + (kx) + & (ee eka 4+ ce bd tnt, 
oY 31 a 

1 


1 1. 
Cy ae a a «» = 2°71B28 2%. 
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Trigonometrical Series 


: 1 ae : 
sin x =X oe tee +. ..ad inf. 
eax el oe ad inf. 
7 a Mga: : 
tanx =x+4x°+4+ + Bt} Sea. . ed ink: 
15 315 


Maclaurin’s Theorem 
1 


f&® =f@ + xf'@ +5 5 xf'O +. sete — ix" (0) +. 


Taylor’s Theorem 
Single Independent Variable 


SOE +N) =f) MPC) + FHF") +. FG HSOCD 
ie aha 
Two Independent Variables 


fethy +h =feoy + (h Z +) 
+o (SE + ome ZL +HoL) 4, 
oy (oy 


9.9. MAXIMUM AND MINIMUM VALUES 

By Taylor’s Theorem, neccesary conditions for a maximum or 

a minimum of f(x, y) are 2f, ws = of/ey = 0. 

Provided that —4 ots ed > (<= f)" then f(x, y) is a maximum 
ax** dy? ~ \axdy/ ’ ; 

or minimum ee as 6°f/0x? and 0?f/éy? are both negative 

or both positive. 
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Section 10 
INTEGRAL CALCULUS 


10.1. STANDARD INTEGRALS 


ff) Sf dx 

x x"*1/(n + 1) except when n = —1 
1/x logex 

1/(ax + 6) * log. (ax + b) 
ekz 1 ekz 

at? a*/k log. a 

sin (ax + b) —_ * cos (ax + b) 
cos (ax + b) ; sin (ax + b) 
tan (ax + b) “ log sec (ax + 6) 
sec? x tan x 

cosec? x — cot x 

sec x tan x sec x 

cosec x cotx —cosecx 
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f[@ Sf dx 


1/\/(a? — x?) — sin=* (x/a) 
1/VG2 + 4%) log - [x + V(x? + a] = sinh-1(x/a) 


VG? = a2) log= [x + V2 — @)] = cosh-1(x/a) 


1/(x? + a?) : tan~1(x/a) 


aa) log (= —£) (> a) 
of 
1/(a? — x?) — el (2 i i 4 (x <a) 
2a a-x 
sinh x cosh x 
cosh x sinh x 
sech? x tanh x 


f'OMFOD log f(x) 
OMS) 2 fd) 


Integration by Parts 

fuS ax = uv — fotar 

Hence f log (x) dx = x log x — x 

10.2. REDUCTION FORMULAE 

J cos" 6.40 = + sin 0 cosr- i¢4¢% nat cost 26 d0 


| sinv 6 ao = — Xcos Osint=*0 +” + “A =t/ sin"~? 6 d@ 
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J sin™ 6 cos" 6 dé 
—1 


L = sin” * O.cost*3 6 + = = J sinn- 2 Ocos" 0 dé 


sin™*+ @ cos"-1 @ ae 4 _ sin” 0 cos"~? 0 dé 
m+n 


1 
m+n 
| tan 0 a0 = J 

ni— 


i tan*-1 06 — | tan"~? 6 d0 


| cot" 6 a0 = fs ! cot? 6 — f cot"? 0 do 


10.3. DEFINITE INTEGRALS 
271 /p 27n/p : 
li sin pt dt = I cos ptdt =0 ninteger 
1 /2p 
a sin ptdt = [ cos pt dt = 1/p 
ai sin? pt dt = eae cos? pt dt = 7/4p 


hs cost ade = X=! 


C cos"~? 6 dé 


pr sin" 0 d0 = td sin"-? 6 0 


/2 er 
™ sin” 0 cos" 0.d9 = ™@—! [7 sinm-2g cos",0 d0 
m+n 
- /2 
nah Sol Ng sin” 0 cos"~*,0.d0 
m+n Jo 


If p and q are integers, p # q, 
27 . 27 , 
I; sin pt sin gt dt = (} sin pt cos gt dt 
Qa 
= ip cos pt cos gt dt = 
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If p = q these become respectively. 
c sin? pt dt = 7 Ay 4 sin 2pt dt = 0 


{* cos? pt dt=7 
0 


10.4. AREA BOUNDED BY A CURVE; LENGTH OF 


CURVE 
ds y= f(x) 
Cartesian Coordinates Ss 
Area A = |" f(xde 
1 
A si ) "3 /( + pdx 
copie Z v( te 6) xy dx Xo ad 
Polar Coordinates 62) 


Area A = iy jr? dO 
41 
pera | in Vir? + (dr/d6)?]d0 
eA 
3 1" VU + 7? (d6/dr)?Jdr 
1 
10.5. VOLUMES AND SURFACES OF REVOLUTION 


When revolved about Ox, the curve y =f) generates a 
volume and a surface of revolution. 


Volume = ip my*dx 
zy 


Area = [2m ds = [2 VU + yadx 
7 aid x 


10.6. SIMPSON’S RULE FOR AREA AND VOLUME 


Area Under a Curve Let the area be divided up into an even 
number n of equally wide strips (width h) by ordinates y,, 
Jas+ + + Yn+1 then 
A = $h[(\1 + yn+r) + 20s + ys t+ yr t+...) 

+ 4(y2 + Ya Jess )) 
Volume of a Solid Let the solid be divided up into an even 
number n of slices by (n + 1) parallel planes at equal distances 
h apart and let Aj, Ao,. . . Anyi be the areas of the sections of 
the solid then 


V = th [(Ay + Ans) + 2 (Ag + As + Az Se 
+4(4,.+44+46+.. .)] 
10.7. CENTRE OF GRAVITY 


The coordinates (<,¥) of the centre of gravity of an area are 


od Ne ses 
ia oy & j= 


10.8. MOMENTS OF INERTIA 


Let I = Mk? be the moment of inertia of a body of mass M 
about a given axis; then k, which has the dimension of a length, 
is called its radius of gyration about that axis. 

Let Jz, Joy be the moments of inertia of the plane area A of 
uniform density o about Ox, Oy, then 


Z2 z2 
los ™ I oty® dx Toy = i ox*y dx 
T1 Ty 


Theorem of Perpendicular Axes (applicable to plane laminae) 
Let Ox, Oy be a pair of perpendicular axes lying in the plane of 
an area A and let Oz be an axis perpendicular to both, then 


Toz = Toz + Toy 
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Theorem of Parallel Axes Let OO’, GG’ be a pair of parallel 
axes at a distance h apart and let GG’ pass through the centre 
of gravity of a body of mass M, then 


Too: = Teg: + Mh? 


Let the solid of revolution and the surface of revolution have 
uniform densities p and o respectively then 


Solid Ios =p |  Anytdx Toy =p i my*(ty? + x2)dx 
Ty 1 
Surface Toz => of 2my%ds =e of” ary? /(1 oe y)dx 
11 
Toy = of” 2ry(dy? + x*)V/(1 + y)dx 
1 


10.9. THEOREMS OF PAPPUS AND GULDIN 


(i) Arc AB of length s does not intersect Ox. 
G is centre of gravity of the arc. 
Surface generated by the revolution of AB about 
Ox is s X 277. 


(ii) Area A does not intersect Ox. G is the centre 
of gravity of A. 
Volume generated by the revolution of A about Ox 
is A X 2nj. 
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Section 11 
DIFFERENTIAL EQUATIONS 


11.1. FIRST ORDER EQUATIONS 


Equation Solution 


(a) Variables Separable 


oo Bae Pd 
fe) 2 = $0) rem 
d : 
fy) 4 = (x) Ir (y)dy = J $(x),dx + C 
Examples 
(i) Decay Function: 
dx|dt = — kx xm Con* 
(ii) Tension in a belt just 
slipping over a pulley: 
dT/d0 = wT T = T, e#8 
(iii) Resisted motion under g 


gravity. Resistance propor- v = k (1 — e7*4) 
tional to velocity v. 


dv/dt = g — kv 


eee 


Equation Solution 


(iv) Current 7 in a circuit (L,R), 
t seconds after a voltage Eo 
is applied to the circuit. 


3s + Ri= Ep, jaa — em FtlLy 
dt R 


(6) Solution by Integrating y = exp (—JPdx) x 
Factor [fQ exp (fPdx)dx + C] 


oe +Py=Q Integrating factor exp ({Pdx) 


(P, Q are functions of x (This method may be applied 
only) to (iv) above) 


(c) Homogeneous Equations 
sD 
P > +Q=0 


(P, Q homogeneous 

functions of x,y) 

Put y = vx so that 

dy/dx = v + x dv[dx 

and equation takes the Solve as (a) above 


form x2 = ¢(v) 
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11.2. SECOND ORDER EQUATIONS 


Equation 


(a) d*y/dx? = f(x) 


(6) d®y/dx? = f(y) 


itp = 9, 


Example 


Inverse Square Law 


dx?/dt? = — yu/x? 


given v = 0 when x =r, 


t=0. 


(c) Simple Harmonic Motion 


Multiply by integrating 


factor 2dx/dt 
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Solution 


y =Sif() dx + Ax +B 


dy|dx)? 


p= 
= WO) dy + A 


©. Sdy]V (2ffQ) dy + A] 
=+x+8 


v = dx/dt $0? = u(- ie 
we 


1 1)28 
dxjdt = VQu)(+— +) 
& F 


Now solve as 11.1 (a) 


(dx/dt)? = n?(a? — x?) 


(dx/dt = 0 when x = a) 
dx/dt = + n(a? — x?)!/? 
x = acos (nt + «&) 


@) 


Equation Solution 


Examples 
(i) pr Pendulum: nae 


S > +& sin 6= = A cos {[V(g/D] ee s} 


e as +¢ 6 = 0 (approx.) Period T = 27v/(i/g) 


if 0 is small. 


(ii) Compound Pendulum: A 
and k are respectively the 
distance of the centre of 
gravity from the axisand = @ = 4 cos {[/(gh/k?)]t + e} 
the radius of gyration 
about the axis. Period T = 271/(k?/gh) 


7 _ sin? = 0 
ee 0 
© 

if 0 is ole 


+35 gh 5 4 = 0 (approx.) 


(d) % + 20% + by =0 


dx 
(a, b constants) 
Put y = Ae** 


A? + 2aA+b=0 
(auxiliary eqn.) 
Ax, Ag are —a + +/(a? — b) Continued on p. 56 
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Equation Solution 


Case I. a? >b y =e-*{A exp [V/(a? — 5)x] 
A, 42, real and unequal. + Bexp [—V/(a@? — 6)x]} 
If b =0, y = A + Be~?% 
Ifa=0,b6 = —-—k <0, 
y = Roeosh kx + S sinh kx 


Case II. a? = b. 

Ax, 42, real and equal. y=(R+Sxe" 
Case III.a* — b = —p? <0. 

A1, ha, complex. y = Re~% sin (px + é) 


Conjugates —a + ip. 


11.3. LINEAR DIFFERENTIAL EQUATIONS OF ORDER 
n WITH CONSTANT COEFFICIENTS 


which may be written 
(an.D" + an-1 D®-*+...+a4D+a)y=0orf(D)y=0 (i) 


where D = d/dx. 
Put y = Ae**, Auxiliary equation: f(A) = 0 which determines 
n values of A so that 


n 
y= >» A,e*+? 
=1 


r 
4 
Example: Whirling of shafts. os —nty=0 


Auxiliary equation: /* — n* =0 
y = Ae™ + Be~™ + Reosnx + S sin nx 
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11.4. PARTICULAR INTEGRALS 
The general solution of 
f(D)y = $@) (ii) 


is y = u + v where u is the general solution of (i) the comple- 
mentary function and v is a particular integral of (ii) given by 


wine 
v = FD) $(x) 


y = Complementary function + Particular integral 


Particular integrals may be found by the use of the following 
operational results. 


ae az 1 C2 ae. Ge 
f(D)e = [@e FD e 7@ e 
f(D) | (ax + b) = f(—a) SM} (ax +8) 


1. on 1 ‘9 
f(D?) cos f(—a?) Cos 


f(Dyle*V(x)] = e% f(D + a) V(x) 


(ax +6) = | (ax +6) 


11.5. FURTHER SECOND ORDER EQUATIONS 


Equation Solution 

a*s a ‘ > 
@ 7 tx=e x = asin (nt + €) + c/n 

(c = const.) 
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Equation Solution 


() & + nex = asin pt x = asin (nt + e) 
(p # n) 
P + =a pt 
The case of resonance: xX = asin (nt + €) — od cosnt 
p=n 2n 
d?x : 
(c) — + 2a he bx Complementary function 
dt? given by 11.2 (d). 
= c sin pt Particular Integral is 
c sin (pt — «) 


V[(6 — p?)? + 4a’p?] 
where tan ¢ = 2ap/(b — p?) 


(d) Homogeneous Equations 
oD tare 2 Pls d(x) 


Ms x = ef es 


ee 
dx at 
ady* _d’y _ dy 
a di? at 
d? 
a t@-Ds ® + by 
= ¢(e') Solution as in 11.4. 
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Section 12 
STATISTICS 


12.1. DEFINITIONS 


Statistical Data A set of measured (numerical) values of x 
called the variable or variate (continuous or discrete). 


Population The totality of all individual values of the variable. 


Sample A finite set of observed values drawn from the popu- 
lation. 


Frequency Distribution Let the range of the variable be divided 
into n class intervals. Let N be the total frequency of individual 
values and let n, be the frequency in the ith class, then 


n n 
N=2Xn, _ Relative frequency f, = x ufi=1 
i=1 i=1 


Statistical data may be displayed by (a) frequency and relative 
frequency curves, (6) histograms, (c) frequency polygons. 


Mode Value of x for which the frequency is a maximum. 


Median Value of x at which the population is divided into 
two numerically equal parts. 


Quartiles Three values of x which divide the population into 
four equal groups. 


Percentiles Values of x which divide the population into 100 
equal groups. 
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12.2. SAMPLE STATISTICS 


Random sample of size N. Number of class intervals n. 


n n 
rth moment about the origin = m,’ = < a ax’ = DY fixe 
i=1 i= 
Mean X = 1st moment about the origin = m,’ = & fix 
l 
rth moment about the mean = m, = W Dn, (x, — X)" 


= Xfi (x, ~ #)" 
m, = 0 and mz = s? = Xf, (x, — ¥)? = U fix? — 
m2 = Sample variance s = Sample standard deviation 
12.3. POPULATION PARAMETERS 


Similar measures for the whole population are 


#4 = population mean o? = population variance 
o = standard deviation. 


For a continuous variable x having a relative frequency or a 
probability density function (p.d.f) y = f(x), 


+o 
rth moment about the origin = yu,’ = | as x’ f (x)dx 
+0 
Mean = Ist moment about the origin = 4,’ = p = | sae xf (x)dx 


+ 
rth moment about the mean = yy, = yen (x — mw)" f(x)dx 
Population variance = fig = 0? = fla’ — pu? 
+o 
=|"? @— sonar 


12.4. ELEMENTARY PROBABILITY THEORY 


An event comprises an outcome or a group of outcomes of an 
experiment. Outcomes may be represented by sample points in 
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appropriate coordinates. Sample space comprises the set of all 
such sample points. 
Let p; @=1, 2, . . ., m) be the probabilities assigned to n 
sample points, then & p; = 1. 
i 


Probability of Event A = P(A) = Xp, 
A 


the sum being taken over the probabilities of the set A of those 
sample points associated with the occurrence of the event A. 


Compound Probability Gi 
P(A U B) = probability of the occurrence of the event A or B 
or both. 


P(A © B) = probability of the simultaneous occurrence of 
events A and B. 


P(A U B) = P(A) + P(B) — P(A 2 B) 


Events A and B are mutually exclusive if the occurrence of one 
precludes that of the other, i.e. A © B = ¢. Therefore 


P(A U B) = P(A) + P(B) 


If Ai, Ao, As, . . . are mutually exclusive events, the probabi- 


lity of the occurrence of at least one is & P(A)). 
i 


Conditional Probability 


The probability of the occurrence of event B subject to the 
condition that event A has occurred is a conditional probability, 
denoted by P(B/A). 


P(B/A) = P(A 1 B)/P(A)_ or 
P(A CX B) = P(A) P(B/A) = P(B) P(A/B) 


If events A and B are independent, P(B/A) = P(B) so that 
P(A A B) = P(A) P(B) 
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If the events Ai, Az, Az, . . . are independent, the probability 
of their simultaneous occurrence is 


P(A) P(A2) P(Az). . 2 = II P(A,) 


Bayes’ Theorem 

Suppose an event B can arise from any one of the initial or 
“causal” events Ai, Ao, . . ., An. The probability that a par- 
ticular causal event A, has resulted in the event B is given by 


P(A,|B) = P(A,) P(BIA,) / PIG ) P(BIA) 


f(x,y) is the Joint Probability Density Function of x and y 
if f(x,y) is the probability of the simultaneous occurrence of 
values x and y. 


PLS rnaee 
eo Jaw LY) ax dy = 1 
If x and y are independent variables 
fy) = fied fy) 
Product Moments of a Bivariate Distribution 
With respect to x = y = 0, 
: +o f+ 
ta’ =f .. xty" f (x, ydx dy 
With respect to the means , and my of x and y, 
+0 [to 
te =f [*? = wo" — mad'flesyax dy 
If x and y are dependent (correlated) variables 


Covariance = cov (x,y) = fii # 0. 
Correlation coefficient = p(x, y) = M11/0,0, 


6,, 6, being the standard deviations of x and y. 
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12.5. VARIANCE OF A LINEAR COMBINATION OF 
RANDOM VARIABLES 


If z = La,x,, the mean Z = pw, = Day. 
i i 
var Z = a,” = Lao? + D1 La,4yp 4400; 
i ce | 


if the variables x, are correlated. 
= Xa,o/? if the variables x, are independent. 
i 


12.6. DISTRIBUTION OF SAMPLE MEANS 
Random sample, size n — x,(@i= 1, 2,. . ., n) 
Sample mean = < = : Lx; 
i 
Mean of sample means = wz = 
Variance of sample means = (0%)? = 0?/n 
12.7. THEORETICAL FREQUENCY DISTRIBUTIONS 
(a) Binomial Distribution (Discrete Variable) 
Probability of x “successes” in n trials = f(x) = (") rg 
where p = probability of “success’ and g = (1 — p). 
Mean = np. Variance = npq. 


q+py"= 3 (7) pgs = 3 fo) 
z=0 z=0 


In N samples, each of 7 trials, the frequencies of 0, 1,2,... 
n successes are terms of the expansion of N (q + p)". 


The proportion of successes in n trials = z = x/n. 


p.d.f. of z = (fA pregna-2 


>’ 
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Mean of z = * (np) = p. 


Variance of z = 5 ("pq) = pq|n. 


(b) Poisson Distribution (Discrete Variable) 


Derived as an approximation of the Binomial Distribution. 
Let m = np be the mean number of successes in 7 trials and let n 
become large and p = m/n small. 


Probability of x successes = e~™ m*/x! 
Mean = Variance = m 


(c) Normal Distribution (Continuous Variable) 
Denoted by N[u, o], Mean yu. St. dev. o. 


we —( — »)? 
pdf. = fs) = — = exp {=e 


(-0 <x< +o) 

Standard Normal Variable = u = (x — y)/o 

Normal Probability Integral for the calculation of cumulative 
normal probabilities 


1 i u me : 
Van . exp (—}u") du is tabulated. 
12.8. ELEMENTARY SAMPLING 


(a) If samples are drawn from a normal population, the sample 
means are normally distributed. 

(b) Let x; be N[v1, 01] and x2 be N[u2, o2]—a pair of inde- 
pendent variables—then a,x; + 2X2 is 


N[aipr + deofle, / (41201? + ag?027)] 
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(c) Central Limit Theorem Sample size n drawn from a popu- 
lation of mean yp, st. dev. o. 


Sample mean ¥ and y = (% — p)(/n)/o 
p.d.f. of y > N[0, 1] as n— o0 


(d) Confidence Limits for the mean y of distribution N[pu, oJ 
from a sample, size n, mean x. 


FA1960//n (95%) #42.580//n (99%) 


12.9. SMALL SAMPLE DISTRIBUTIONS 


(a) Small sample y,(i = 1, 2, . . ., n) drawn from population 
N[{O, 1). 


n 
u= 4° = X y? has a x? distribution f(y?) with n degrees of 
i=l 
freedom (d.0.f.) given by 


Bi be (y?)t"~ te~ 4x2 
S02) = fu) = FS 


The cumulative x? distribution i fu) du has been tabulated. 


(6) Confidence Limits for the variance o? of a normal variable. 
Sample, size n, variance s*. (100 — «)% confidence limits for 
o” are 


ns?/¥,? > 02 > ns?/yq? 
where P(z? < 41) and P(y? > x2?) are each 30%. 


(c) F-distribution If u, v are independently distributed y? vari- 
ables with 7, vg d.o.f. then if F = u/v, + v/ve, the p-d.f. p(F) is 


F 
known and the cumulative distribution function ‘ P(F)dF has 


been tabulated for various d.o.f. 


65 


(d) F-ratio Test 


To test, on the basis of variance, whether two samples, sizes 1, 
Ng, Variances 5S”, 52”, are likely to have been drawn from the 
same population of variance 0”. 


u = n4S,2/07, v = NeS2?/o” have x? distributions with 
(ny a 1), (ng — 1) d.o.f. 


mS,? N2S2” ire 
= = 6,7/6.? 
a -tl &—) ee 


where 6,2, 62? are unbiased estimates of o?, on the null hypothesis 
that the samples are drawn from the same population. 


(e) t-distribution 
If u is N[O, 1] and v has a x? distribution with » d.o.f., the 
probability distribution of t = uv/»/+/v has been tabulated. 
Sample size n drawn from N[p, o]. 
Taking u = (*% — b)V/n[o [see 12.8 (c)] and v =ns?/o? 
having a x? distribution with n — 1 d.o.f. 
t= (% — wV(n — 1)/s has a ¢ distribution with n — 1 d.o.f. 


(f) Confidence Limits of the Mean 


- S ms s Aer O 
> Seat ee wae [1001 —« %] 


where P(t > tx)2) is «/2. 


(g) t-test: Difference of Two Sample Means: Is there any essen- 
tial difference between two similar samples of data? Suppose 
x,y are distributed N[., o,] and N[u,, o,) respectively. Two 
random samples, sizes n;,”,, means X,/. 

On the assumption that the samples are drawn from the same 
population, x and y are distributed with common variance 
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o = 0, = 0, [justify the assumption by use of F-test, 12.9(d)] 
and fl, = My. 


=~ EP tf NyMy ) ya 
t= 5 poesia 4 mar with (n, + ny, — 2) d.o.f. 


and 6 is the unbiased estimate of o from pooled data. 


12.10. GOODNESS OF FIT 


Population divided into k mutually exclusive classes. f,, m, are 
respectively the observed and expected frequencies in class i. 


k ead 2 
ee has oy distribption with (k — 1d 
jung 


12.11 LINEAR REGRESSION 


Bivariate distribution. Sample size N. Sample points: (x;, y,). 
Frequency f, and N = & f; 


(a) Least Squares Regression Line y = a + bx 
- ee apy | 
a =f — bx where t = wuSixed = wsin 


NU fix — UfixDfyi 


sis eilnstias 5 fat — Opn K 


(b) Unbiased Estimate of Population Correlation Ratio p 


is given by 
r=s' |s 's/ = Zfix: — %) (yi — J) 
ww ABs — ZR BA — IVP 
where 5,’ Sy’, Sry’ are sample st.devs. and covariance respectively. 


(c) Significance of r On the hypothesis that p = 0, 
r/[(n — 2)/G — r?)] has a f¢-distribution with (n — 2) d.o.f. 
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If the sample value of t is significant, the hypothesis is rejected 
and a degree of correlation is accepted. 


(d) Rank Correlation (Spearman) Items ranked and assigned 
ordinal numbers in respect of some characteristic (usually 
magnitude). 

Values of x have ranks u, and values of y have ranks v. 


62d,” 


Spearman’s coefficient r= 1 — 
i (n° — n) 


where d, = u; — v4. 
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Section 13 
NUMERICAL ANALYSIS 


13.1. ROUNDING OFF ERRORS 


Let A, be the rounded-off measured value of a quantity whose 
exact value is A,. 
(Maximum rounding error = } unit of the last decimal 
place in A,) 
Absolute error = a = A, — A,. 
Absolute error modulus of A, = |a}. 
Relative error = a/A, 


Let a1, a2, ag be the absolute errors of three measured values 
(A,)1, (Az)a, (Ay)3 respectively, then if 


(A,)s = (A,): of (A,)2 |as| < |a,| 7 |a2| 


(A ra my (A rt = (A re a3 Pe 
or (A,)1 X (A,)2 (A Mi (Ay)s (A,)2 

A,)2 = {(Ay i} oa 

(A,)2 = {(A,)1} Fe IPl - a re 


13.2. ERRORS IN FUNCTIONS 


(a) If u = f(x) and x is subject to error, then 


bu = ZL 5x 
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where du, 6x are the absolute errors in u, x respectively. 
(6) Ifu = f(x, y,z.. .) then 


=ZoxrtZoyt+ Zor. 


13.3. TABULAR VALUES OF f(x); DIFFERENCE 


NOTATION 
Constant tabular interval = h 
Value of x Value of f(x) Notation 
xX-2 fo — 2h) f-2 
X-1 fo —/A) Fe 
Xo (x o) f o 
X1 f (x ° + h) f 1 


X2 fo + 2h) fr 


Forward Difference Operator A 

First difference: Afn = Af (xn) = f(%Xn + h) — f(x) 
Second difference: A(Afn) = A?fn = Afny1 — Afn 
rth difference: A'fn = A'~*faar — Af 


13.4. LINEAR INTERPOLATION 
If x, = Xo + ph, then f, = fo + p Afo 
13.5. SHIFT OPERATOR E 


E{ fo} = f@ + A’) E*{ f)} = ELE fOOH 
= E{ f(x + h)} = f(& + 2h) 
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In general 

E°{ f(x)} = f(x + ph) (p +ve, —ve or fractional) 
A more concise notation for this is 

Ef, = fna+p 
In particular, Ef-1 = fo .. ET Fg fos 

E*? fy = fija = f(x + 4h) 


13.6. RELATIONSHIP BETWEEN A AND E 
Afn =fn+1 — fn oe (1 + A)fn = fn+1 = Efn 
“ E=1+A and £’=(1 + Ay 
(r +ve integer) 
13.7. GREGORY-NEWTON FORWARD DIFFERENCE 
INTERPOLATION FORMULA 
Value of f(x) required at x = x9 + ph. 
ho = E°fo = (1 + A)?fo 


—1 —1\(p-—2 
={f + pAfo +Pe- ) arg, PP we ) Asp, 
tr 


13.8. BACKWARD DIFFERENCE OPERATOR V 
Defined by Vf(Qin) = fn) — f On — A) or Vn = fn — fr-2 
and Vif, = Vf, — VH3fq-1 

E~*fn = fa-i = fn — Vin = (1 — V0 fn 

E-+=(1—V) or E=(1—V)7? 
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13.9. GREGORY-NEWTON BACKWARD DIFFERENCE 
INTERPOLATION FORMULA 

fo = Efe =(1—V*fo= {fot pV fo +249 vy, 
+Pet Pt ve +...| 


13.10. CENTRAL DIFFERENCE OPERATOR 5 


6 is defined by S5fn+112 =fn+1 — fn for all integral values of 
n. The operator 6 may be applied to any column of the difference 
table so that 


a = Ofn+r1/2 ea Ofn-112 

O8fn+1/2 — O*fns1 = Oo ete. 
In particular, 67+? f,;2 = 6*f; —6"fp (r =1,2,3,...) 
Ofn+112 = fn+1 —fr ss EF Fi 41/2 — E~*!?fa 4132 

= (£1? — E-*!?) fn+1j2 

Therefore 6 = E1/? — E-1/2 

E§=(E—-1)=A 6%=A2/E= A211 + A) 
13.11. CENTRAL DIFFERENCE INTERPOLATION 

FORMULAE 


These formulae use information placed symmetrically with 
respect to the interpolation interval. 


(i) Bessel f, = 
fot p ofii2+ ve ) (0%fo + Of) +e ea Ofiia 


+ Pt EP (84f + fi) +. De 
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(ii) Everett 
Using (i) of 13.10, the odd order differences may be eliminated 
from the Bessel formula above giving 


fo = (1 — p)ifo + pf + PP ~ NEP) 524, 
ee ee 


[Further coefficients will be found on p. 56 of Interpolation and 
Allied Tables (1.A.T.), H.M.S.O.] 


(iii) Stirling 
fo = fo + Sil6f—-1)2 + Ofti2) + S26*fo + Ss (6%f-1)2 + O%fij2) 
ie ae 


The coefficients S,, So, . . . will be found on p. 58 of 1.A.T. 


13.12. LAGRANGIAN INTERPOLATION 


f(x) takes the values fo, fi, . . .. fn for x = Xo, X1, ~~. ., Xn 
respectively. It is required to find a polynomial which takes 
these (n + 1) functional values for the same values of x. The 
required interpolate is found by calculating the value of this 
polynomial for the given value of x. 


Lagrangian Interpolation Coefficients 


L(x) = 
Ce AoE — Ri) wi Me om Fie Riss). . A(X — Xo 
(x; = Xo)Mx; — X1) + yee (x; on Xy-1)(X;—- X41) es (x; baad Xn) 


defines the (m + 1) coefficients for i = 0, 1,2,...n. 
L(x) are polynomials of degree n. 
L(x;) = OF xf) and L(x,) = 1. 
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The Lagrangian Polynomial is defined by 
n 
~ LQ) = XL,)f, which is of degree < n 
i=0 


Since L(x,) = f, for i=0, 1, 2,. . ., m, L(x) is the required 
polynomial. If f(x) is of degree <n, it may be shown that 
f(x) = LO). If f(x) =e, 


Ln=c= > Linx k >> L(x) = 1 (check) 
i=0 i= 0 


13.13. AITKEN’S METHOD FOR INTERPOLATION 
WITH UNEQUAL INTERVALS 


f(x) has values fo, fi, fj at x = Xo, x1, Xy respectively (x1 — Xo 
# X;— X1, in general). Determine the interpolate at x. 


(x ios Xo) fo 
(x — %)) fi 
(xy — Xo) 


(x—%*1) for) 
(x — x) fos(x) 


Second approximation = fo1; (x) = Gi x 


First approximation = fo,(x) = 


13.14. NUMERICAL INTEGRATION 
Hf ie ng ph [2 feode =n [fy dp 
T) 


(i) Forward Difference Integration Formula (derived from 
Gregory—Newton formula) 


[P*reoae = h| nf += A fo + 


+ Oo Ae +. ; | 


moar 
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and, in particular, 

[P rooax =h E aah = “5 A +A +.. | 
Simpson's Rule (an even number n equal subdivisions of width /) 
[2 foode = Vl fo + fa + fat fot. + faa) 

+ 4(f. + fot. <). + fe-1)} 
Estimate of error = — = {Atfo + Ati t+... + Atfp_a} 
Trapezoidal Rule (n equal subdivisions of width h) 

[PP foods = th fe + fi + fats. A fad + fa 

Estimate of error = 

— 5 (W8fn 08 fi) + 2 U%fu— 155fo) 


where y is the averaging operator defined by 


Mf> = 4 (fp+1/2 + fp-1)2), 
HOF, = 40"fp+1)2 + Ofp-1;2) (I.A.T. p. 54) 


(ii) Central Difference Integration Formulae 
(a) Derived from Bessel’s interpolation formula. 


i fodde = h {fo + fd — 5g O%fo + Hf) 
a coe ae , 
=h|nfsa- 2 5 Phas + a a5 MOMfuia + . | 
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(6) Derived from Stirling’s interpolation formula. 
ca Mh lice. 4 ) 
[P fends = 2h (fo + Lop a Offot +. 
13.15. NUMERICAL DIFFERENTIATION 
dp _1 
If = - 4 = 4 he sa 4 
X=Xotph f(x) f'(p) zf P) 


(i) Central Difference Formulae 
(a) Derived from Bessel’s Interpolation Formula. 


£0) =1 |Sfarat ACp — IN6*fo + fr) + Gp? — 3p + Dar 
+55 Op — DG" =p = NOYfo + 5471) +. «| 
At half-way points f’(%o + 3h) 
=F lite A fin to Hint. J 
(6) Derived from Stirling’s Interpolation Formula. 
FG) =F [Af-210 + fara) + po%fo + 4 Gp? — 1N6°F-11 
+ fay) + 5 pp? — I) d*fo +. . 
At a tabular point xo, p = 0 so that 
f00) =4 [NOf-ai0 + fire) — 5 O10 + Hf) 
+ a OF -a0 + Of) +. . ‘ 
= 7 {udfy — 2 pdf +d wdfet...| 
h 6 30 
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(ii) Forward and Backward Difference Formulae (derived from 
the Gregory-Newton Formulae) 


(a) Forward Difference formula. 


(0) = 4 |Afo + AOp — DA%o + 23p* — 6p + DAY, 


+98 


For the tabular point xo, put p = 0 so that the Forward Differ- 
ence formula becomes 


fo) =; {af — A%, + 4A% — 2th + 1 Ate soil 4 


(6) For a point near the end of the table, use the Backward 
Difference formula 


fo) =" MZ + 4V%fy + AV + 1V%fy +3V% ro: y 
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